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I. 


INTRODUCTION 


Previous  work  on  the  problem  of  diffraction  by  a  nlane  angular 
sector  (see  Fiqure  11  has  been  with  the  exact  eiqenfunction  solution 
[1,?]  but.  results  have  been  presented  for  only  the  auarter  olane  M.e., 
k2  =  1/2).  The  eigenfunction  solution,  of  course,  requires  a 
larqe  number  of  terms  to  obtain  a  converqed  result  and  this,  as  well 
as  the  complexity  of  the  solution,  restricts  its  useful Iness.  As 
stated  by  Satterwhite  [?],  "the  solution  is  qeneral  and  needs  somp 
simol if ication  before  it  can  hecome  a  useful  tool  for  diffraction 
studies." 


The  ourpose  of  this  paner  is  to  make  a  simol if ication  in  the 
exact  eigenfunction  solution  so  that  the  solution  may  be  more  easily 
used.  Not  only  are  the  eigenvalues  found  for  the  quarter  plane  as 
has  been  done  previously,  but  the  first  few  contributing  eigenvalues 


1 


are  found  for  almost  any  anqular  sector.  These  results  are  presented 
in  qraohical  form  for  easy  use. 

The  second  purpose  of  this  paper  is  to  combine  the  eigenfunction 
solution  with  UTD  wedqe  diffraction  theory  to  obtain  the  current  on 
the  anqular  sector  usinq  only  two  or  three  eiqenfunctions.  The  use 
of  only  a  few  eiqenfunctions  is  possible  because  the  eiqenfunction 
solution  is  used  to  find  the  current  only  in  the  reqion  very  close 
(i.e.,  reO.lx)  to  the  vertex.  The  UTD  is  then  used  to  determine  the 
current  everywhere  else.  By  contrast,  Satterwhite  [?]  used  fifty 
eiqenfunctions  to  calculate  the  current  out  to  only  one  wavelength 
from  the  vertex. 

Results  will  be  shown  for  current  distributions  on  angular  sectors 

p 

of  various  shapes  (i.e.,  0.1<k  <0.91  and  results  will  also  be  shown 
for  the  far  field  when  the  anqular  sector  is  illuminated  by  a  short 
dipole.  The  angular  sector  will  then  he  used  alonq  with  superposition 
and  the  UTD  to  find  the  current  on  four-sided  thin  plates. 

II.  THE  SPHERO-CONAL  COORDINATE  SYSTEM 

The  sphero-conal  coordinate  system  is  one  of  the  few  coordinate 
svstems  in  which  the  vector  wave  equation  is  separable  and  it  is  in 
this  svstom  that  the  plane  anqular  sector  is  a  coordinate  surface. 

The  sphero-conal  system  is  described  by  three  coordinate  surfaces 
which  are  a  sphere,  an  elliptic  cone  and  an  elliptic  half-cone  as 
illustrated  in  Figure  2.  A  plane  angular  sector  is  the  e=n  coordinate 
surface  of  the  sphero-conal  system  and  its  angle  is  determined  by  the 

p 

ellipticity  parameter  k  of  the  system.  The  r,n,<j>  coordinates  are 
related  with  the  cartesian  x,y,z  by  the  following  equations: 

2  p~ 

x  =  rcosejl  -  k1  cos  <j> 


? 


v  =  rsin0sin<t> 


z  =  rcos^1 


fis 


? 

COS  0 


m 


w Here  k'?  =  l-k?,  0  <  k?  £  1,  0  <  0  <tt  and  0  <  4>  £  ?v .  The  0  and  ^ 
anoles  are  not  the  same  as  the  anqles  of  the  spheriea1  coordinate 
system.  A  ooint  on  the  angular  sector  (Figure  ))  has  coordinates 
rs,es,<t>s  where  0$  =  tt  and  <t>s  =  cos"^  (-^p--).  The  element  ds  of  the 
surface  on  the  angular  sector  is  given  by: 


ds  =  rsdr$ 


k1  sin<)>s 

J  1  -k '  ?cos^ft>s 


d<t>s 


The  surface  0=m  and  <b=0  as  well  as  the  surface  0=ir  and  d>=ir  result 
in  half  planes.  For  a  more  detailed  discussion  of  the  sphero-conal 
system  see  Satterwhite  [?  ]  and  Kraus  and  Levine  [T,4]. 


II! .  EXALT  SOLUTION  FOR  FIELDS  AND  THE 
ANGULAR  SECTOR  CURRENT  DISTRIBUTION 

The  T  field  due  to  the  current  density  U  is  determined  by  the 
dyadic.  Green's  function  which  we  must  determine.  The  T  field  expression 
is  of  the  form 

?(R)  =  jujfi  f  r(R,R').J(R'ldv  .  (31 

v 

H  is  a  field  point,  IT'  is  a  source  point  and  v  is  the  source  region. 

In  Satterwhite' s  paper  [1]  we  can  see  that  the  Green's  function  is 
given  by 


r(R.R')  = 


j<  T. 

g 


^2  (R)  "Jx  <*') 

V 


!a[ 


(TT* ) 


ql 


f«) 


for  [IT]  >  0?'] 


t 


and 


rfR.IT')  = 


J  <  l 
9 


^ 

V 


N1,  fRl  N1!  f R * ) 
ql  ;  q! 

Aql 

for  [IT]  <  [IT1] 


(5) 


where:  <  is  the  wave  number 

g  is  the  symbolic  index  which  means  the  class  of 
the  odd  or  even  solution  of  the  problem 

A  is  the  normalization  constant. 
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The  vector  wave  functions  are 


JK 


WKr’ 


9? 


sin?fl+k  ,?sin?<t> 


.[/ l-k?cos?9  - 


7l-k ,?cos?<t>  0  (*) 


„  Z  t  fieri 

v-v  v,,jaW  V6' ♦>)?* 

[rZw^i(<r)l’  / - 9 - O- 

+  ~  vgl -  -  [>/l-!< 1  cos  <t>  0nlfft)  *Md>)  S  + 

r?  ?  ?  ?  9  - 

r^k^sin  0+k'^sm  * 

+  Jl- k2cos?q  O’, f A )  *  ,(d>)p]  .  f 71 
ql  q  I 

The  Z v (icr )  function  is  the  spherical  Bessel  Function  where  the  super¬ 
script  II  means  that 

7vf.crl  =  h^ficr)  ffll 


and  the  suDerscript  I  means  that 

V<r)  =  Jv(<r)  .  (9) 

The  0(a)  and  <j>(<b)  functions  are  Lame  functions  defined  by  the  scalar 
wave  equations 

J l-k2cos2A  ^  (J  l-k2cos2A  +  [(v(v+l)k2sin2e+y]e  =  0  ( in i 

and 

/1-k,2cos?*  ^  ( Jl-k’?cos^  ^1  +  [  (  +1)k,?sin?'b-u>  =  n.  Ml) 

The  two  Equations  MOl  and  Ml)  are  coupled  throuqh  the  eiaenvalues 
v  and  y.  The  H  field  is  qiven  by  the  curl  of  the  T  field  as  follows: 

W1  =  L.  vxf(ff)  =  -v  x  f  T(K,ff'  1.3(R‘  )dv  (121 

u>y  -V 

and  the  current  of  the  anqular  sector  is 

■J(JT)  =  hxfrf(r,A=TT,0«t><iO  -  H(r,fi=|TT,TT«{><?TT)  ]  .  (131 

IV.  EIGENVALUES  AND  THE  0(A)  AND  <S>(i>)  FUNCTIONS 

From  the  above  expressions  we  can  see  that  we  always  need  the 
first  even  and  odd  0(A)  and  4>(<b)  functions  and  the  eigenvalues  associ¬ 
ated  with  them.  Rv  using  the  infinite  continued  fraction  method  as 
in  the  first  three  pairs  of  the  main  contributing  eigenvalues 

u  and  v  were  found.  These  are  plotted  in  Figure  3  where  the  ordinate 
is  k  and  thn  ahscissa  is  either  (v  or  y)  for  0<v<1.5. 

The  qeneral  expressions  for  the  fl(Pl  and  <P(<p)  functions  are 
Vfil  =  z  Am  cos(?m  -  ^)A  fm=. n, M41 


6 


0Q(  ft)  =  r.  Am  sin (2m  -  a 

fm--. . .  ,-2,-1 ,0, 1 ,2, . . .  1 

M6i 

<J>ef  (t>)  =  T  Bm  cos  m* 

fm=0,l ,2,3, . . . ) 

(16) 

<t>Q( a)  -  r  R  sin  m<b 

(m=l ,2,7, . . . )  . 

(17' 

Knowinq  the  eigenvalues  u  and  v,  the  coefficients  A„  and  R  mav  he 
_  mm 

found  for  0.1<_k  <0.d  by  the  continued  fraction  method  used  in  [  1,?] 

—  -  ? 

for  the  auarter  plane  (k  =  1/2)  case.  The  procedure  is  tedious  and 

laborious.  The  results  are  presented  in  Fiqures  A_Q.  The  coefficients 

A  and  R  mav  be  determined  from  the  fiqures  with  sufficient  accuracy 
mm 

to  use  in  Fquations  (12)  and  MR).  Thus,  with  the  aid  of  the  fiaures, 
the  current  in  a  reoion  close  to  the  vertex  mav  rather  easily  he 
determined  for  almost  any  anqular  sector  as  we  shall  see  in  the  followino 
sections. 


V.  E  FAR  ZONE  FIELD  APPROXIMATION  FOR  UNIT 
DIPOLE  SOURCE  LOCATED  CLOSE  TO  THE  TIP 


An  investigation  of  the  eigenvalues  v  and  u  shows  some  inter- 

p 

estinq  results.  First  we  note  that  for  any  k 


Vlel  <  Vlo2  <  V?el 


These  three  eioenvalues  v-j^,  v-|0p  and  v,,  ^  have  values  which  varv 

from  0  to  1.6.  If  we  examine  Equations  (61  and  (7),  we  see  that  the 

[r7  (■cr)]1 

-  for  <r<<l  correspondina  to  the  wave 


<r 


functions  7  (<rl  and 
functions  m  and  become 

/ir^Krlv 


7vfKrl  =  .iv(icr)  = 


?v+^  rrv+1.5) 


MR  1 


R 
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Figure  9.  The  B  coefficients  of  the  ?(<j>)  function  for 
the  second  even  solution. 


ItZ^kt)]'  [rjv(tcr)]' 


2v(?v+l  )rfv+0.5) 


For  small  values  of  v  the  expressions  HP)  and  (IP)  are  relatively 
larger  than  they  are  for  hig  values  of  y.  This  means  that  in  the 
dvadic  Green's  function  the  first  few  terms,  say  three,  have  relatively 
hiq  values  comoared  to  the  others.  So,  when  we  want,  to  qet  an  approxi¬ 
mate  expression  for  the  far  zone  electric  field  in  the  case  of  a  unit 
dinole  source  located  close  to  the  vertex,  we  can  use  only  the  first 
few  eiaenvalues  in  the  expression  of  the  dvadic  Green's  function  and 
obtain  accurate  results.  In  Figure  10  we  see  two  locations  of  the  unit 
dipole  source.  In  Figures  11  to  15  we  see  the  far  zone  electric  field 

O 

with  0.1<k  -0.9,  corresponding  to  Fig.  10a.  In  Figures  16  to  20  we  see 
the  far  zone  electric  field  for  the  second  location  of  the  unit  source 
dipole.  The  computations  in  Figures  16  to  20  are  for  angular  sectors 
with  0.15k^50.5. 


Figure  10, 


e  Ee  far  zone  field 
ose  to  the  tip  of  p 
ource  on  x  axis  and 


Fiau^e  1?. . 


The  %  far  zone  field  In  the  x=o  plane  due  to  a 
unit  source  dipole  close  to  the  tip  of  plane 
angular  sectors  for  k^=0,l  to  0,9  (source  on 
x  axis  and  parallel  to  y). 
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90° 


0  0.2  0.4  0.6  0.8 


Figure  13.  The  E$  far  zone  field  in  the  x=o  plane  due  to 
a  unit  source  dipole  close  to  the  tip  of 
plane  angular  sectors  for  k2=0,l  to  0,9 
(source  on  x  axis  and  parallel  to  y). 
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SOUR 


far  zone  field  on  the  plane  angular  sectors  due  to  a  unit  sojrce 
close  to  the  tip  for  k*=0,l  to  0,9  (source  on  x  a\i-s  ard 


VI.  CURRENT  OISTR  HUTION  ON  THE  ANGULAR  SECTOR 


From  Equation  MSI  it  can  he  shown  that  the  current  distribution 


TWH  =  -  3cY„  cJ-— -  'RJ.T'Rjdv.ft  ,  (-7I 

0  i  V  '  07  0  0  0? 


,i  (<r)  [rj  fKr)]'  ?  ? 

E-V V”  *o><*'*  +  — f" - TW" 


Jv  (rr) 


for  r  <  r 
-  o 


Tf(tf)  =  .  2vcYo  T  ;  wj?(^0.3(^0)dv  o0?m 


h(?\<r)  [rh^Vtcrl]'  9 

r  t  vo?  .  / ,  no  .  07  1-k'  cos  a  ,  . 

[_Vfvo?+1)  - icr - +  - - !TrTm* -  VUVj+  a 


b(J]()cr) 

*7!  -fkr- 


Q 1  f  »  *  1 

TT5T-  «eT*'rj 


for  r  >  rQ  . 


The  expressions  indicate  that  the  current  in  the  r  direction  has  the 
form  l/sin<J).  This  means  that  the  current  on  the  edges  qoes  to  infinity. 
But  we  know  the  current  distribution  due  to  geometrical  optics  and 
the  edge  diffracted  field  does  not  give  infinite  current  on  the  edaes 
when  the  incident  electric  field  is  normal  to  the  edge.  So,  we  conclude 
that  the  l/sin*  must  arise  because  of  the  tip  diffracted  field. 
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Near  the  tip  the  current  distribution  is  larqely  determined 

by  the  spherical  Ressel  function  .1  <rl  when  the  source  is  within 

v 

a  distance  rQ  >  r.  For  <r«l  the  expression  of  the  current  will  depend 
on  the  j^(icr)  which  have  small  values  of  u.  As  mentioned  earlier, 
the  small  eiqenvalues  are  the  first  two  or  three.  These  qive  the 
dominant  current  near  the  tip.  With  the  help  of  Equation  f13)  for 
the  even  eigenfunctions  the  current  distribution  is 


^el  =  8(<r) 


iVel  !el_ 


and  for  the  odd  eigenvalues  it  is 


^o?  =  A(<r) 


vo?_1 


/  .  1  Jl-k  ,?cos?<b  . ,  ,  v- 

V0?  k'sin<j>  *o?.^r 


where 


R  =  -  ?kY. 


^el  (v) _ 

V  i  +  l 

?  r(vel+i.suel 


f  N^fR’l.^'ldu 


A  =  -  2kY  - ?Z—  - - — -  r  M^fR’l.'JfR'ldu 

’  0?  r  (vo2+  1.5)  «o2 


Yq  is  the  free  space  admittance.  From  the  expressions  of  the  current 
near  the  tip  we  can  see  in  general  that  the  current  depends  on  three 
different  functions. 


In  Figures  IF,  17  and  1R  we  can  see  the  form  of  those  functions 

p> 

for  an  angular  sector  with  k=0.3.  Far  from  the  tip  6c r> 1 )  the  tip 
diffracted  current  must  be  given  by: 


Jtd(R) 


-jicrr  <l>-j  f  d») 


r  + 
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Figure  16a.  The  J  =(kr) 


Figure  16b.  The  Jf  current  distribution  as  in  Fig.  16d.  as  seen  from  a  different  aspect. 


Figure  18a.  The 


Figure  18b.  The  J  current  distribution  as  in  Fig.  18a  as  seen  from  a  different  aspect. 


wHe>-e  the  functions  ♦  are: 


=  A0+A|COS<t>+Aj,cos/>,4'+. . .  07 ' 

=  B0-*-B i S i i 


T  e  unknown  coefficients  A-  and  8^  can  be  found  bv  usina  the  EFIE 
to  generate  two  systems  of  linear  equations  using  point-matchinc.  '"-e 
EFIE  is 

n  x  fnC  (r)  =  n  x  /  ( -a>?pe3sG+7;  G }dS 

where  G  is  the  free  space  Green's  function  and  in  this  problem 
is  given  by 


7s  "  ^total 


7s  ■  *»  *  *  ?T0 


0  _<  <r<l  f  see  Eo.  ( ?01^j 

V 

<r  >  1 


in' 


or  in  words  for  .r.l,  the  current  is  the  sum  of  the  geometrical 
optics  current,  the  diffracted  current  from  the  edges  and  the  tio 
diffracted  current  respecti vely. 

VII.  ANGULAR  SECTOR  FAR  FIELD  RESULTS 

By  using  the  EFIE  with  just  three  sample  points,  we  found  the 
three  unknowns  in  the  tip  diffracted  current  expressions  in  Equations 
(27)  and  (28).  The  integration  of  EFIE  was  terminated  20  wave¬ 
lengths  from  the  sample  point  since  the  integrated  functions  go 
quickly  to  zero  with  distance.  The  sample  points  are  typically 
located  about  10*  from  the  tip.  In  Figures  19,  20,  and  21  we  can  see 
the  total  field  for  an  infinitesimal  dipole  located  near  the  angular 
sector  and  normal  to  it.  The  distance  of  the  dipole  from  the  tip 
is  i cr  =  2.0  with  0Q  =ir  and  <J>0  =  */2. 


_ _ . ,  -  —  -  -•  - 
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•'in.  FINITE  THIN  PLATES 

A  finite  thin  plate  has  ho  n  edges  and  vertices.  We  can  use 
the  angular  sector  solution  in  section  VI  to  obtain  the  currents  in 
the  regions  close  to  the  vertices  and  the  UTD  to  obtain  the  currents 
elsewhere.  If  the  vertices  are  far  apart  by  more  than  a  wavelength 
or  so,  we  can  neglect  the  interaction  between  the  vertices  and  use 
the  expressions  in  Sections  VI  and  VII  plus  superposition  to  obtain 
the  total  contribution  from  the  vertices.  That  is,  we  use  the  EFIE 
as  in  the  previous  section,  except  that  the  number  of  sample  or  match 
points  will  be  n  times  the  number  of  unknowns  used  in  the  anqular 
sector  case. 

Results  for  a  "U  hy  TX  souare  pla*°  shown  in  Fiqurp  ??. 

Hi  •>  to  the  nature  of  our  solution  wp  spe  a  larqo  "spike"  of  current 
at  thn  four  corners  in  Fioure  9,a.  The  current  in  Fiourp  ??a  mav 
Hp  comnared  with  that  Drpspn‘pr)  hv  Ko  anr1  *'’ffra  [  f]  ,  and  it  will 
hp  seen  that  t^e  results  arp  siqmficantly  different,  at  the  four  carners 
since  their  results  show  no  singularity  of  the  current  at  the  corners 
for  this  size  plate.  We  believe,  that  because  the  interaction  between 
corners  is  small  (or  negligible)  for  large  plates,  it  is  inevitable 
that  a  strong  singularity  be  observed  in  current  there.  This  con¬ 
clusion  is  the  opposite  of  that  in  [5].  It  may  be  that  the  solution  in 
[5]  is  not  converged  for  large  plates. 

Figure  23  shows  the  current  on  a  3X  by  3x  plate  that  is  not  square, 
but  in  the  shape  of  a  rhombic  with  two  corners  characterized  by 
k^=0.3  and  the  other  two  by  k^=0.7.  In  both  Figs.  22  and  23  the  edge 
condition  of  the  current  is  not  shown  for  purposes  of  clarity  nor  is  it 
used  in  the  far  field  computation  in  the  next  paragraph. 

Finally,  Figure  24  shows  the  far  field  pattern  of  a  short  monopole 
or  stub  radiator  at  the  center  of  a  8x  square  ground  plane.  The  result 
is  compared  with  an  unpublished  result  by  Marhefka  who  obtained  his  re¬ 
sult  using  the  UTD  with  vertex  diffraction.  The  two  results  are  seen  to 
be  in  very  close  agreement. 


30 


Figure  19.  Total  electronic  field  for  five  different 
angular  sectors. 
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Figure  20.  Total  electric  field  for  five  different 
angular  sectors. 


Figure  21.  Total  electric  field  for  five  different 
angular  sectors. 
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IX.  SUMMARY 


In  this  paper  we  have  presented  In  easy  to  use  graphical  form 
the  first  few  contributing  eigenvalues  for  any  angular  sector  having 
a  wide  range  of  included  angle.  We  have  shown  that  it  is  only  these 
first  few  eigenvalues  that  contribute  to  the  current  in  a  region  very 
close  to  the  tip  itself.  Beyond  this  region  the  GTD  may  be  used  to 
obtain  the  current  thereby  overcoming  the  need  to  use  large  numbers 
of  eigenvalues  and  eigenfunctions. 

We  then  applied  the  combined  eigenfunction-GTD  technique  to 
find  the  currents  and  resulting  radiation  fields.  To  demonstrate 
the  validity  of  the  solution,  the  radiation  pattern  of  a  monopole 
at  the  center  of  a  square  plate  is  compared  with  that  calculated  with 
an  unpublished  GTD  solution  which  Includes  vertex  diffraction.  The 
two  independent  results  are  in  very  close  agreement. 


short  monopole  in  the  center. 
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APPENDIX 


The  following  tables  give  the  lower  order  eigenvalues  for  angular 
vectors  with  0.1  <  <  0.9  and  for  0.0  <  v  <  1.0  and  0.5  <  v'  <  1.5. 
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